Deconfinement and criticality in extended two-dimensional dimer models 
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A square-lattice hard-core dimer model with links extending beyond nearest- neighbors is studied 
using a directed-loop Monte Carlo method. An arbitrarily small fraction of next-nearest-neighbor 
dimers is found to cause deconfinement, whereas a critical state with r~^ distance dependence 
of the dimer-dimer correlations persists in the presence of longer dimers preserving the bipartite 
graph structure. However, the critical confinement exponent governing the correlation of two test 
monomers is non-universal. Implications for resonating-valence-bond states are discussed. 
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Dimer models have a long history in classical statistical 
physics More recently, they have also emerged 

as central models in modern theories of strongly corre- 
lated quantum matter, e.g., high-temperature cuprate su- 
perconductors and frustrated antiferromagnets 0, 0j 
1^ JMi JJjJ ■ The dimers then represent singlet-forming 
electron pairs. In order to model the quantum fluctua- 
tions of the dimers and realize a short-range version of 
Anderson's resonating valence bond (RVB) state ,12.. l-^ . 
Kivelson, Rokhsar, and Sethna introduced a Hamiltonian 
with a term flipping (resonating) pairs of parallel nearest- 
neighbor dimers on the two-dimensional (2D) square lat- 
tice 13 . The purely classical dimer model retains it rele- 
vance also here: It was shown that the equal- weight sum 
over all dimer configurations is the ground state of the 
Hamiltonian when the resonance strength —k equals the 
potential energy cost v of each resonating pair of dimers 
(the RK point) 0. This state is critical; the dimer- 
dimer correlations decay with distance as and two 
inserted test monomers are correlated with each other 
as r~^/^ Q. As it turned out, away from the RK point 
the dimers form long-range order and the monomers are 
exponentially confined 0, Q, [l5'|. Hence this system 
does not give rise to the desired RVB state with no bro- 
ken lattice symmetries and deconfined monomers (corre- 
sponding to spin-charge separation 13]). Moessncr and 
Sondhi recently showed that a true extended RVB phase 
does appear in the quantum dimer model on the trian- 
gular lattice 13 • Following this insight, a large body 
of work has been carried out in order to characterize 
classical an d q uantum dimer models on various lattices 
H US El El mill- Moreover, there are currently 
intense activities in gauge theories related to quantum 
dimer models H H Efl O El El . 

To date, research on dimer models has focused mainly 
on planar lattices, i.e., ones that have no intersecting 
links. This class of models can be analytically solved 
(in the form of Pfaffians) with the aid of a theorem 
by Kastelyn j^, and thus the quantum ground states 
at the corresponding RK points are characterized as 
well. It has often been stated that the inclusion of 
valence bonds (dimers) extending further than between 



nearest-neighbor sites will not change the physics pro- 
vided that the probability of longer bonds decreases suf- 
ficiently rapidly |6||,|17,]. However, the fact that there are 
qualitative differences between bipartite and nonbipar- 
tite lattices, e.g., the square and triangular cases men- 
tioned above, does raise the question of potentially im- 
portant effects of short bonds connecting two sites on the 
same sublattice of the square lattice. Such bonds will in- 
evitably appear in realistic systems away from the limit- 
ing cases m 1^ represented by the nearest-neighbor 
(iVi) dimer models. Introducing next-nearest-neighbor 
(iV2) links along one of the diagonals makes the square 
lattice equivalent to the triangular one. It has already 
been shown that an arbitrarily small fraction of such di- 
agonal dimers destroys the critical iVi square-lattice state 
and leads to deconfinement [HI ITsj l , as in the isotropic 
triangular lattice The model with links along both 
diagonal directions, i.e., the full 2D square lattice with 
iVi and A''2 bonds, is not solvable by Kastelyn's theorem 
Intuitively, one might suspect that the critical state 
is immediately destroyed in this case as well, but no cal- 
culations have been carried out thus far. One might also 
speculate that introducing longer bonds between the two 
sublattices does not lead to deconfinement, but there are 
no results available to back this up. 

In this Letter, two extended dimer models on the 
square lattice are studied — the nonbipartite model with 
A'^i and N2 dimers as well as the bipartite lattice with 
iVi and A^4 (fourth-nearest-neighbors, of which there are 
eight per site). An efficient directed-loop jl^ Monte 
Carlo algorithm is used to sample the full space of hard- 
core dimer configurations, with fugacities Wi assigned to 
the different dimer types. The results confirm that a 
low, most likely infinitesimal, concentration of N2 dimers 
leads to deconfinement, and that the presence of A'4 
dimers does not. However, while the dominant dimer cor- 
relations are ~l/r^ also when W4 > 0, those correlations 
involve N4 dimers; the A^i correlations decay as a higher 
power of l/r. The structure in wave- vector space also 
becomes much more complicated. Furthermore, the crit- 
ical confinement exponent governing the monomer cor- 
relation function is nonuniversal, changing continuously 



from —1/2 in the pure Ni case to —1/9 (to high numer- 
ical accuracy) for the pure N4 model. 

Directed-loop algorithm. — This algorithm is an adap- 
tation of a quantum Monte Carlo method with the 
same name, in which updates of the system degrees of 
freedom are carried out along a self-intersecting path 
at the endpoints of which there are defects not allowed 
in the configuration space contributing to the partition 
function. When the two defects meet they annihilate, 
the loop closes, and a new allowed configuration is cre- 
ated. The conditions for detailed balance in the process 
of stochastically moving one of the defects are expressed 
as a coupled set of directed-loop equations [2^. The ap- 
plicability of this scheme to dimer models was first re- 
alized by Adams and Chandrasekharan [2^ (an alterna- 
tive cluster algorithm was constructed in Ref. Here 
an algorithm for the multi-length dimer problem will be 
presented; a simplifying representation of the dimer con- 
figurations will also be introduced. The algorithm will 
be described only for the case of the N1-N2 model, but 
the scheme applies directly to any range of the links. 

The links connected to a given site are labeled as shown 
in Fig. n^a). For a dimer configuration on a periodic 
Lx L lattice, each site is numbered according to the type 
of dimer it is connected to. For any given site, the other 
member of the same dimer can then be easily found. The 
central object in the directed-loop algorithm is a vertex, 
which in this case consists of a site and all its surrounding 
sites to which it can be coupled by a dimer. The state 
of the vertex is the number (1-8) assigned to the central 
site. A step in the directed-loop algorithm is illustrated 
in Fig. ^b). The vertex is entered through the dimer, 
and one of the seven other surrounding sites is chosen as 
the exit. The dimer is then flipped from the entrance to 
the exit, the central vertex site remaining connected to 
it. The exit site already has another dimer connected to 
it, and this site now becomes the entrance in the next 
step of the algorithm. In each step a defect is hence 
moved one link ahead, leaving behind a healed dimer 
state. The first entrance site is chosen at random and the 
corresponding defect remains stationary. It is annihilated 
when the moving defect reaches it, whence a new allowed 
dimer configuration has been generated. In the present 
case the defects are monomers, and the intermediate two- 
monomer configurations can thus be used to measure the 
correlations between two inserted test monomers. 

The key to an efficient algorithm of this kind is that 
the probabilities for random selection of the seven possi- 
ble exit sites can be chosen in such a way that detailed 
balance is satisfied without any further accept/reject cri- 
terion. Each step then moves one dimer, and a full loop 
can accomplish very significant changes to the dimer con- 
figuration. The directed-loop equations [i^ give the con- 
ditions for detailed balance in terms of weights ajk for 
the processes in which a vertex in state j is entered at 
site j and exited at k (transforming the vertex into state 
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FIG. 1: (a) Labeling of the links of the N1-N2 dimer model, 
(b) A step in the directed-loop update, in which the entrance 
to the vertex is at site 3 and the exit is at site 8. The thick 
bond indicates the location of the dimer. 



k). The actual probabilities Pjk — ajk/wj, which implies 



Detailed balance is satisfied if a 



jk — ^kj • 



In the N1-N2 model, the vertices can be classified as even 
(e) or odd (o) according to the numbering of Fig. ^a); 
there are then four weights: Oee, Qoo, Oeo, Q-oe- In princi- 
ple, one can include "bounce" processes where j = k, but 
in the present case they can be excluded. Including only 
the seven no-bounce exits, the directed-loop equations 
reduce to 



Wi 



4a„ 



W2 = Sflee + 4aeo, 



(la) 
(lb) 
(Ic) 



This system is underdetermined and has an infinite num- 
ber of postive-definite solutions. Here the following solu- 
tion will be used: For wi >W2, 



while for 1x12 >wi, 



Oo 



doe = a-eo = 'i«2/7, 

[wi - 4u;2/7)/3, 



floe = fleo = Wl/7, 

{W2 - 4u;i/7)/3. 



(2a) 
(2b) 



(3a) 
(3b) 



There is no guarantee that this is the best solution, but 
the resulting algorithm performs very well and allowed 
for studies of lattices with '^10^ dimers. The fugacity of 
the A^i-bonds is set to unity (except in the case of the 
pure N2- and iV4-models) . The proCTam was tested using 
known results for the pure Ni case |j| and by comparing 
with local Metropolis simulations for small lattices. 

Results. — Dimer-dimer correlations in the full close- 
packed system and monomer-monomer correlations in 
the system with two test monomers will be discussed. 
The monomer correlations M(r) were obtained by ac- 
cumulating the distances between the stationary and the 
moving monomer in the directed-loop update. As has be- 
come customary |23, the normalization M(r = 1) = 1. 
Several types of dimer-dimer correlations can be defined. 
Here Dy, (r) will be defined in the following way: If site i 
is connected to a dimer in the set E, a variable s{i) = 1, 
otherwise s(i) — 0. The correlation function is then 
D-^{vij) = {s{i)s{j)). Results will be presented for cases 
where E contains a single dimer or half of the dimers of 
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FIG. 2: Dimer (a) and monomer (b) correlations along the 
direction (r, 0) for the Ni-N2 model at W2 — e"*, compared 
with those of the pure A^i model. The solid lines show the 
agreement with the known power-laws for the A'^i model. 
Note that for the pure A'^i model M(r) = for even r. The 
results were obtained using L = 1024 lattices. 



one type. For the N1-N2 model, the correlation functions 
Di, D2, Da, and Db are thus defined corresponding to 
the sets {!}, {2}, A = {1, 3}, and B = {2, 4}. Analogous 
definitions are used for correlations Z?i, of A^i dimers 
and D4, Djj of A'4 dimers in the N1-N4 model. 

In Fig.|2Ia), dimer correlations for the nonbipartite 
N1-N2 model with W2 = e"** (corresponding to a concen- 
tration p2 ~ 0.3% of N2 dimers) are compared with those 
of the pure A'^i model. There is a clear deviation from the 
decay, showing that the N1-N2 model is not critical 
at this very low concentration of dimers. As shown 
in Fig.[2Ib), the test monomer correlation approaches a 
nonzero constant, i.e., the system is deconfined in con- 
trast to the critically confined A^i model. The very signif- 
icant changes seen already at a very low concentration of 
N2 dimers suggest that an arbitrarily small concentration 
indeed causes deconfinement. 

Turning now to the bipartite N1-N4 model, its dimer 
correlations I?^(r) (TVi dimers) and D£,{r) {N4 dimers) 
are compared with DA{r) of the pure Ni model in Fig. 01 
Surprisingly, the two correlation functions decay with dif- 
ferent power-laws: Dd{t) apparently always follows the 
same form as Da^j) of the pure iVi model, whereas 
DA{r) decays faster once W4 > 0. In Fig. 0] the Fourier 
transforms of the correlation functions involving a single 
type of dimer are shown for the pure A'^i and N4 mod- 
els. For the A^i model, the dominant correlations are at 
q = (tt, 0). A logarithmic divergence of the peak height 
with the system size can easily be observed (not shown 
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FIG. 3: Various dimer correlations along the direction (r, 0) 
in the Ni-N^ model (on L = 1024 lattices). The fugacities 
W4 = e~* and e~^ correspond to concentrations p4 ~ 0.065 
and 0.17, respectively. The solid line shows the asymptotic 
form r~^. 



here). The N4 model exhibits more complicated corre- 
lations, with very broad peaks that show almost no size 
dependence up to the largest size (L — 2048) that was 
studied. Nevertheless, a logarithmic divergence should 
eventually occur here as well, considering the real- 
space correlations. The details of the nature of these 
critical correlations remain to be elucidated. In the A^i- 
A'4 model, the wave-vector structures shown in Fig. ^ 
change with W4. Most notably, the Di peak at (tt, 0) is 
suppressed as the fraction of N4 dimers increases, and a 
nondivergent structure at (tt, tt) emerges. 

An unchanged dominant critical dimer exponent might 
suggest that the monomer correlations should also remain 
of the pure iVi form r^^^^ . This is not the case, however. 
Results for M(r) at r = L/2 — 1 are shown multiplied 
by L" in Fig. |5fa). Here a is adjusted to give a flat L 
dependence, and hence the critical form M{r) ~ r~" is 
extracted. The known a — 1/2 is used for the pure A^i 
model. For the pure N4 model the exponent is consistent 
with a — 1/9 (to an accuracy of 1%). To show that the 
exponent changes from 1/2 already at a low concentra- 
tion of N4 bonds, results for W4 = e~^ {p4 w 0.9%) are 
scaled with a = 1/2 in Fig. Efb). This scaling clearly 
fails, and is instead consistent with a « 0.485. 

Conclusions. — The results obtained here demonstrate 
that the 2D square-lattice dimer model becomes decon- 
fined when a very low concentration (likely infinitesi- 
mal) of non-bipartite (next-nearest-neighbor) dimers are 
introduced. The system remains critically confined in 
the presence of longer bipartite (fourth-nearest-neighbor) 
dimers, but the corresponding two-monomer exponent is 
nonuniversal. In contrast, the dimer-dimer correlation 
exponent appears to remain unchanged, although the na- 
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FIG. 4: Fourier transform of the dimer-dimer correlation Di 
of the Ni model (left) and D4, of the A'4 model (right) calcu- 
lated on L = 128 lattices. 



ture of the dominant correlations changes. 

These findings are relevant to quantum dimer mod- 
els as well. Various resonance terms can be introduced, 
and corresponding RK points can then be demonstrated 
in the same way as has been done for other models 
(^OiEHl' Hence, it is clear that an RVB state can be 
realized on the square lattice once next-nearest-neighbor 
dimers are allowed. On the other hand, the emergence of 
new structure in the dominant dimer-dimer correlations 
on the bipartite graph including fourth-nearest-neighbor 
links indicates that phase transitions between different 
long-range ordered valence-bond-solid states can be real- 
ized. Such quantum order-order transitions have recently 
been discussed by Vishwanath, Balents, and Senthil [23 ■ 
Extended quantum dimer models on the square lattice 
should thus have very rich phase diagrams and may also 
be relevant in the context of the cuprates. 

This work was supported by the Academy of Finland, 
Project No. 26175. The calculations were carried out 
on the Condor system at the University of Wisconsin- 
Madison. 
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FIG. 5: Finite-size scaling of monomer correlations in the A'^i- 
model, (a) The best scaling for three different cases, (b) 
Failure of a scaling with the pure A^i exponent a = 1/2 for 
the N1-N4 model at low fugacity w^. 
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